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Abstract 

Recent breakthroughs in cosmology reveal that a quarter of the Universe is composed of dark 
matter, but the microscopic identity of dark matter remains a deep mystery. I review recent 
progress in resolving this puzzle, focusing on two well-motivated classes of dark matter candidates: 
WIMPs and superWIMPs. These possibilities have similar motivations: they exist in the same 
well-motivated particle physics models, the observed dark matter relic density emerges naturally, 
and dark matter particles have mass around 100 GeV, the energy scale identified as interesting 
over 70 years ago by Fermi. At the same time, they have widely varying implications for direct and 
indirect dark matter searches, particle colliders, Big Bang nucleosynthesis, the cosmic microwave 
background, and halo profiles and structure formation. If WIMPs or superWIMPs are a signifi- 
cant component of dark matter, we will soon be entering a golden era in which dark matter will 
be studied through diverse probes at the interface of particle physics, astroparticle physics, and 
cosmology. I outline a program of dark matter studies for each of these scenarios, and discuss the 
prospects for identifying dark matter in the coming years. 

PACS numbers: 95.35.+d, 13.85.-t, 12.60.Jv, 04.65. +e 



1 




I. INTRODUCTION 



In recent years, there has been tremendous progress in understanding the Universe on the 
largest scales. Observations of supernovae, the cosmic microwave background (CMB), and 
galaxy clusters have provided three stringent constraints on and ^a, the energy densities 
of matter and dark energy in units of the critical density. These results are consistent and 
favor (f2M,^A) ~ (0.3,0.7), as shown in Fig. ^ The amount of matter in the form of 
baryons is also constrained, both by the CMB and by the observed abundances of light 
elements together with the theory of Big Bang nucleosynthesis (BBN). Although there are 
at present possibly significant disagreements within the BBN data, the CMB and BBN 
data taken as a whole are also impressively consistent, providing yet another success for the 
current standard model of cosmology. 

Through these and many other observations, the total energy densities of non-baryonic 
dark matter, baryons, and dark energy are constrained to be 0, |3|] 

fioM = 23% ±4% (1) 
Qb = 4% ±0.4% (2) 
^A = 73% ± 4% . (3) 

These results are remarkable. At least two of the constraints of Fig. must be wrong 
to change the conclusions about the central values of flu and ^a. These results are also 
remarkably precise — the fractional uncertainties on all three are (9(10%). Given that just 
a decade ago the range 0.2 < fl-DU ~ 0.6 was allowed and ^a = was often assumed, 
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this represents spectacular progress. Although much of cosmology remains imprecise, as 
we will see, the quantum leap in precision in these three quantities already has dramatic 
implications for particle physics. 

At the same time, recent progress in cosmology is probably best viewed as the first 
steps on the road to understanding the Universe. Consider an historical precedent: in 200 
B.C., Eratosthenes determined the size of the Earth. On a day when the Sun was directly 
overhead in Syene, Eratosthenes sent a graduate student to measure the lengths of shadows 
in Alexandria. He was then able to extrapolate from the known distance between these two 
cities to determine the circumference of the Earth. His answer was 



27ri?e = 250, 000 stadia . (4) 

This result is remarkable. At the time of publication, it was bigger than many expected, 
leading many to be skeptical and helping to earn Eratosthenes the nickname "Beta" J4|. His 
result was also remarkably precise. We now know that it was good to less than 10% |aMl3> 
where the leading source of uncertainty is systematic error from the exact definition of the 
unit "stadion" j8|. At the same time, the achievement of Eratosthenes, though important, 
could hardly be characterized as a complete understanding of the Earth. Rather, it was just 
the beginning of centuries of exploration, which eventually led to the mapping of continents 
and oceans, giving us the picture of the Earth we have today. 

In a similar vein, recent breakthroughs in cosmology answer many questions, but highlight 
even more. Focusing on dark matter, the primary subject of this review, these include 

• What particles form dark matter? 

• Is dark matter composed of one particle species or many? 

• What are dark matter's spin and other quantum numbers? 

• What are its interactions? 

• How and when was it produced? 

• Why does f^DM have the observed value? 

• How is dark matter distributed now? 

• What is its role in structure formation? 

• Is it absolutely stable? 

Although these questions will continue to be sharpened by astrophysical observations 
at large length scales, it is clear that satisfying answers will require fundamental progress 
in our understanding of microphysics. This is nothing new — the history of advances in 
cosmology is to a large extent the story of successful synergy between studies of the Universe 
on the smallest and largest length scales. This interplay is shown in Fig. |2l where several 
milestones in particle physics and cosmology are placed along the cosmological timeline. 
As particle experiments reach smaller length scales and higher energies, they probe times 
closer to the Big Bang. Just as atomic physics is required to interpret the CMB signal from 
t ~ 10^^ s after the Big Bang and nuclear physics is required to extrapolate back to BBN at 
t ~ 1 s, particle physics, and particularly the physics of the weak or Fermi scale, is required 
to understand the era before t ~ 10~® s, the era that contains the answers to many of our 
most basic questions. 
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FIG. 2: Milestones in particle physics and cosmology along the cosmological timeline. 



Here I will review recent progress in dark matter, focusing on new proposals for and 
studies of dark matter at the Fermi scale 

Mf ~ 100 GeV . (5) 

Mf is the scale of the weak interactions; its importance for particle physics was realized by 
Fermi some 70 years ago. Remarkably, it is now appreciated that incisive studies of this scale 
may also yield deep insights into the dark matter problem as well. In the next few years, the 
Tevatron and Large Hadron Collider (LHC) will play the crucial role of opening the door to 
the Fermi scale. These will be followed, we hope, by detailed studies of new physics at these 
colliders and the proposed International Linear Collider (ILC). If dark matter particles do 
in fact have masses at the Fermi scale, these colliders, in conjunction with cosmological and 
astroparticle probes, will be essential for studying dark matter and unveiling its identity. 

I begin with a review of the motivations to consider dark matter candidates with Fermi 
scale masses. These candidates may be divided into two classes: WIMPs and superWIMPs. 
I consider these in turn, discussing the basic scenarios, their implications for particle physics 
and cosmology, and the prospects for identifying dark matter and answering many of the 
other questions listed above. Further background and discussion on both these and other 
dark matter topics may be found in Refs. [9|. 



4 



II. WHY THE FERMI SCALE? 



The particle or particles that make up most of dark matter must be stable, at least 
on cosmological time scales, and non-baryonic, so that they do not disrupt the successes 
of BBN. They must also be cold or warm to properly seed structure formation, and their 
interactions with normal matter must be weak enough to avoid violating current bounds 
from dark matter searches. The stringency of these criteria pale in comparison with the 
unbridled enthusiasm of theorists, who have proposed scores of viable candidates with masses 
and interaction cross sections varying over tens of orders of ma gnit ude. In roughly the 
order in which the y w ere proposed, these include axion s (]H 11. 1^ |. thermally- produced 
gravitinos Q Q El H El El 111 113 , ^utralinos |2ll2j"^axinos 0, Q balls H, 

Kaluza-Klein 



wimpzillas |2q| . self-interacting dark matter [2^ , annihilating dark matter j^iij , i 
dark matterjisl 129^ . branons [3^ 31 1, superWIMPs [H 33j . and many others. 

Candidates with Fermi scale masses have received much of the attention, however. There 
are at least four good reasons for this. First, these proposals are testable. Second, new 
particles at the Fermi scale are independently motivated by attempts to understand elec- 
troweak symmetry breaking. Third, these new particles often "automatically" have all the 
right properties to be dark matter. For example, their stability often follows as a result of 
discrete symmetries that are necessary to make electroweak theories viable, independent of 
cosmology. And fourth, these new particles are naturally produced with the cosmological 
densities required of dark matter. 

The last motivation is particularly tantalizing. Dark matter may be produced in a simple 
and predictive manner as a thermal relic of the Big Bang. The evolution of a thermal 
relic's number density is shown in Fig. |21 In stage (1), the early Universe is dense and hot, 
and all particles are in thermal (chemical) equilibrium. In stage (2), the Universe cools to 
temperatures T below the dark matter particle's mass m^, and the number of dark matter 
particles becomes Boltzmann suppressed, dropping exponentially as e~'^^l'^ . In stage (3), 
the Universe becomes so cold and dilute that the dark matter annihilation rate is too low 
to maintain equilibrium. The dark matter particles then "freeze out," with their number 
asymptotically approaching a constant, their thermal relic density. 

More detailed analysis shows that the thermal relic density is rather insensitive to 
and inversely proportional to the annihilation cross section: ^dm ~ (ca"^)"^- The constant 
of proportionality depends on the details of the microphysics, but we may give a rough 
estimate. On dimensional grounds, the cross section can be written 

GAV = 1 or v') , 6 

mr, 

A 

where v is the relative velocity of the annihilating particles, the factor f ^ is absent or present 
for S- or P-wave annihilation, respectively, and terms higher-order in v have been neglected. 
The constant ai is the hypercharge fine structure constant, and k parameterizes deviations 
from this estimate. 

With this parametrization, given a choice of k, the relic density is determined as a function 
of m^. The results are shown in Fig. |21 The width of the band comes from considering both 
S- and P-wave annihilation, and from letting k vary from | to 2. We see that a particle that 
makes up all of dark matter is predicted to have mass in the range ~ 100 GeV — 1 TeV; 
a particle that makes up 10% of dark matter, still significant (with respect to its impact 
on structure formation, for example), has mass m-^ ~ 30 GeV — 300 GeV. There are 
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FIG. 3: Left: The cosmological evolution of a thermal relic's comoving number density. Right: A 
band of natural values in the {my^, Qy,) plane for a thermal relic IsT 



models in which the effective k is outside our illustrative range. In fact, values of k smaller 
than we have assumed, predicting smaller m^, are not uncommon, as the masses of virtual 
particles in annihilation diagrams can be significantly higher than m^. In any case, the 
general conclusion remains: particles with mass at the Fermi scale naturally have significant 
thermal relic densities. For this reason, a thorough exploration of the Fermi scale is crucial 
if we hope to identify the particle or particles that make up dark matter. Even null results 
from LHC and ILC searches for dark matter are important, as without them, it is unlikely 
that we will be able to exclude the possibility of dark matter at the Fermi scale. 

Given these results, many theories for new physics at the Fermi scale contain promising 
dark matter candidates. The candidates that exploit the tantalizing numerical "coincidence" 
shown in Fig.|21may be grouped into two classes: WIMPs and superWIMPs. In the following 
sections, we consider each of these two cases. 



III. WIMPS 



Weakly-interacting massive particles (WIMPs) have weak-scale masses and weak-scale 
interactions. They are an especially well-motivated class of dark matter particles, and there 
are many examples, including neutralinos in supersymmetry 21 . 2^ |. Kaluza-Klein particles 



in theories with universal extra dimensions 2^1, branons in theories with large extra 



dimensions j30|, |31| , and the lightest T-odd particle in some little Higgs theories [35[ . 

A program of detailed WIMP dark matter studies may be divided into three (overlapping) 
stages: 

1. WIMP Candidate Identification. Is there evidence for WIMPs at colliders from, for 
example, events with missing energy and momentum? Are there signals in dark matter 
search experiments? What are the candidates' masses, spins, and other quantum 
numbers? 

2. WIMP Relic Density Determination. What are the dark matter candidates' predicted 



6 




FIG. 4: Left: Minimal supergravity (mSUGRA) defines a hypersurface in the 105-dimensional pa- 
rameter space of the MSSM. In the studies described here, the underlying parameters are assumed 
to lie on the mSUGRA hypersurface, but deviations in all directions of the MSSM parameter space 
are allowed when evaluating the potential of colliders to constrain parameters. Right: Schematic 
diagram of regions with the right amount of dark matter (shaded) in mSUGRA. This diagram is 
qualitative. The precise locations of the shaded regions depend on suppressed parameters, and 
axis labels are purposely omitted [i^. 



thermal relic densities? Can they be significant components or all of dark matter? 

3. Mapping the WIMP Universe. Combining collider results with results from direct and 
indirect dark matter searches, cosmological observations, and A^-body simulations, 
what can we learn about astrophysical questions, such as structure formation and the 
distribution of dark matter in the Universe? 

Stage 1 is discussed in Ref. fsi^. In the following sections, we will explore how well near 
future experiments may advance Stages 2 and 3. 

To address these issues concretely, it is necessary to focus on one representative example. 
We choose to study neutralinos. Even with this restriction, there are many qualitatively 
different scenarios. A common choice is to study minimal supergravity (mSUGRA), a simple 
model framework that encompasses many different possibilities. In this case, one assumes 
that the underlying supersymmetry parameters realized in nature are those of a point in 
mSUGRA parameter space. In determining the capabilities of experiments, however, it is 
best to relax all mSUGRA assumptions and ask how well the 105 parameters of the general 
Minimal Supersymmetric Standard Model (MSSM) may be determined. This approach is 
illustrated in Fig. HI 

In much of mSUGRA parameter space the neutralino relic density lies above the narrow 
allowed window, and so these possibilities are cosmologically excluded. The regions in which 
the neutralino relic density is not too large, but is still sufficient to be all of dark matter, 
are cosmologically preferred. They have been given names and include the bulk, focus 
point, co-annihilation, and rapid annihilation funnel regions shown in Fig. |3J Results from 
representative models in each of the first two regions are summarized below. For results for 
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FIG. 5: Processes contributing to neutralino WIMP annihilation [isf . 



the other two regions and related studies, see Refs. [221, QSl, Q^, liil • For each model, the 
superpartner spectrum is determined by the computer code ISA JET ji^], and cosmological 
observables, such as the thermal relic density, are determined by the software packages 
DARKSUSY and micrOMEGAs H. 



A. WIMP Relic Density Determination 

To determine the predicted WIMP thermal relic density, one must experimentally con- 
strain all processes contributing significantly to the WIMP pair annihilation cross section. 
This requires detailed knowledge not only of WIMPs and their properties, but also of all 
other particles contributing to their annihilation. This is no small task — the number of 
processes contributing to WIMP annihilation is large, as illustrated in Fig. El All unknown 
parameters at the Fermi scale must be either measured precisely or constrained sufficiently 
so that their effects are known to be irrelevant. Such detailed work relies primarily on 
particle colliders, and we now consider how well colliders may constrain the thermal relic 
density. 
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FIG. 6: Left: The dominant neutralino annihilation process in the bulk region. Right: The 
superpartner spectrum at LCCl, a representative model in the bulk region |4^ . 



1. Bulk Region 

In the bulk region, a much studied model is specified by the mSUGRA parameters of 
Linear Collider Cosmology Model 1 (LCCl): 

LCCl: (mo,Mi/2,v4o,tan/5) = (100 GeV, 250 GeV, -100 GeV, 10) , (7) 

with yU > 0, m3/2 > "^LSP, and nit = 178 GeV. Here mo, M1/2, and Aq are the uni- 
versal scalar, gaugino, and trilinear coupling masses specified at the grand unified scale 
Mqut — 2.4 X 10^^ GeV, respectively, tan (3 = {H^) / {Hf) is the ratio of Higgs boson vac- 
uum expectation values, /i is the supersymmetric Higgs mass, and ^3/2 is the gravitino 
mass. The neutralino thermal relic density at this point is f^^^^ ~ ^-^^ — sig- 
nificantly higher than the range f^^^^ ~ 0.113 ± 0.009 allowed by the latest cosmological 
constraints Nevertheless, the choice of LCCl is convenient, since it has been studied 

in great detail in other studies, where it is also known as Snowmass Points and Slopes Model 
la (SPSla) |46|]. 

In the bulk region, neutralinos annihilate dominantly through XX ~^ ff through a t- 
channel scalar /, as shown in Fig. IHl To achieve the correct relic density, this process 
must be efficient, requiring light sfermions and, since the neutralino must be the lightest 
supersymmetric particle (LSP), light neutralinos. These characteristics are exhibited in 
the superpartner spectrum of LCCl, shown in Fig. IHl It is noteworthy that in this case, 
cosmology provides a strong motivation for light superpartners within the reach of a 500 
GeV ILC. 

To determine the relic density at LCCl, all of the supersymmetry parameters entering 
annihilation processes, including those shown in Fig. [Ul and others, must be determined 
to high accuracy. The LCCl superpartner spectrum makes possible many high precision 
measurements at the LHC. LCCl (SPSla) is in significant respects a "best case scenario" 



9 



6000 I — < — < — < — ■ — , — ■ — ■ — < — < — , — < — < — < — r 




leplon energy E, [GeV] (GeV) 



FIG. 7: Left: The kinematic endpoints of lepton energies from e^e~ l^l~XX provide 

precise determinations of slepton and neutralino masses 47|. Right: Threshold scans may also be 
used to determine slepton masses. In the case of selectron masses, e~e~ threshold scans provide 
higher precision and simultaneously save luminosity [i^ . 



for the LHC The implications of these measurements for cosmology will be summarized 
below. 

The LHC results may be improved at the ILC. For example, superpartner masses may be 
determined with extraordinary precision through kinematic endpoints and threshold scans, 
as shown in Fig. [7| The kinematic endpoints of final state leptons in the process e^e~ — *• 
Z+/~ — > 1^1~XX determine both I and x masses. Slepton masses may also be determined 
through threshold scans. Threshold scans provide even higher precision, and may actually 
save luminosity. This is the case, for example, for selectron mass determinations through 
e~e~ threshold scans, where precisions of tens of MeV may be obtained with 1 to 10 fb~^ 
of integrated luminosity j48fc49ei^J. More generally, the required measurements exploit the 
full arsenal of the ILC, from its variable beam energy, to its polarized beams, to the e~e~ 
option. The results of one study are summarized in Fig. |H| 

The neutralino thermal relic density may be determined by combining the precise deter- 
mination of all relevant supersymmetry parameters and also verifying the insensitivity of 
the relic density to all other parameters. The results depend somewhat on the prescription 
one uses to combine these data. One approach is to choose points in parameter space at 
random, weighting each with a Gaussian distribution for each observable. The relic density 
allowed region is then identified as the symmetric interval around the central value that 
contains 68% of the weighted probability. 

The result of applying this method with 50,000 model parameter points randomly selected 
around LCCl is shown in Fig. El The result is that the ILC may determine the thermal relic 
density to a fractional uncertainty of 

LCCl (preliminary): ^ ^ = 2.2% [^{^xh^) = 0.0042] . (8) 

The current constraint from WMAP and the projected future constraint from the Planck 
satellite are also shown. WMAP and Planck provide no information about the mass of the 
dark matter particle. 
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FIG. 8: The precision with which superpartner masses may be determined at the ILC for charginos 
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masses, the second the constraint from collider studies, and the third the method used to achieve 
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2. Focus Point Region 



In the focus point region, one may choose the representative model, Linear Collider 
Cosmology Model 2, specified by 



LCC2: (mo, Ml/2, Ao, tan/?) = (3280 GeV, 300 GeV, 0, 10) 



(9) 



with /i > 0, 7113/2 > ^Lsp, rrit = 175 GeV. In focus point supersymmetry j51|, |52|) squarks 
and sleptons are very heavy, and so the diagrams that are dominant in the bulk region are 
suppressed.^ Nevertheless, the desired relic density may be achieved [H^, because in the 
focus point region, the neutralino is not a pure Bino, but contains a significant Higgsino 
component. The processes XX ~^ W~^W~, shown in Fig. ^1 and XX ~^ which are 
negligible in the bulk region, therefore become efficient. Neutralino mixing is typically 
achieved when neutralinos and charginos are fairly light and not too split in mass, and so 
the demands of neutralino dark matter motivate supersymmetry with light neutralinos and 
charginos. 

Determination of the thermal relic density in the focus point region requires precise mea- 
surements of neutralino and chargino masses and their mixings. Applying the method de- 
scribed above for converting collider constraints to a constraint on the thermal relic density, 
the thermal relic density may be determined with fractional uncertainty 



LCC2 (preliminary): 



A(^]^/^.2 



2.4% [A{n^h'^) = 0.0026] 



(10) 



"'^ As a result of this property, models like focus point supersymmetry may be challenging for supersymmetry 
discovery and study at the LHC 
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3. What We Learn 

The results of Figs.lQlandliniimply that the ILC will provide a part per mille determination 
of fi-^h^ in these cases, matching WMAP and even the extraordinary precision expected from 
Planck. The many possible implications of such measurements are outlined in the flowchart 
of Fig. HH 

Consistency of the ILC and WMAP/Planck measurements at the part per mille level 
would provide strong evidence that neutralinos are absolutely stable and form all of the 
non-baryonic dark matter. Such a result would at last provide convincing evidence that 
we have produced dark matter at colliders and that we have identified its microphysical 
properties. It would be a landmark success of the particle physics/cosmology connection, 
and would give us confidence in our understanding of the Universe back to neutralino freeze 
out at t ~ 10~^ s, eight orders of magnitude earlier than can currently be claimed. 

On the other hand, inconsistency would lead to a Pandora's box of possibilities, all with 
important implications. If the thermal relic density determined from high energy physics is 
smaller than what is required cosmo logically, these high precision measurements imply that 
neutralinos are at most only one component of cold, non-baryonic dark matter. On the other 
hand, if the thermal relic density determined at colliders is too large, these measurements 
imply that neutralinos must decay (perhaps to superWIMPs — see below), or that the 
neutralino thermal relic density is diluted by entropy production or some other effect after 
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FIG. 10: Left: The dominant neutralino annihilation process in the focus point region. Right: 
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freeze out. 

The implications of LHC precision measurements for the relic density, determined in the 
way discussed above, are also shown in Fig. 121 The LHC precision in the LCCl scenario 
is extraordinary and unusual; for other scenarios, the LHC is unlikely to determine Vly^ 
to better than one or more orders of magnitude. At the same time, even in this "best 
case scenario," the LHC determination of VL-^ leaves open many possibilities. For example, 
comparison of the LHC result with WMAP/Planck cannot differentiate between a Universe 
with only neutralino dark matter and a Universe in which dark matter has two components, 
with neutralinos making up only 80%. Such scenarios are qualitatively distinct, in the sense 
that the possibility of another component with such significant energy density can lead to 
highly varying conclusions about the contents of the Universe and the evolution of structure 
that formed the galaxies we see today. 



B. Mapping the WIMP Universe 

WIMPs may appear not only at colliders, but also in dark matter searches. Direct 
dark matter search experiments look for the recoil of WIMPs scattering off highly shielded 
detectors. Indirect dark matter searches look for the products, such as positrons, gamma 
rays or neutrinos, of WIMPs annihilating nearby, such as in the halo, the galactic center, or 
the core of the Sun. 

If WIMPs are discovered at colliders and their thermal relic densities are determined to 
be cosmologically significant, it is quite likely that they will also be discovered through direct 
and indirect dark matter search experiments. The requirement of the correct relic density 
implies that WIMP annihilation was efficient in the early Universe. This suggests efficient 
annihilation now, corresponding to significant indirect detection rates, and efficient scatter- 
ing now, corresponding to significant direct detection rates. This rough correspondence is 
illustrated in Fig. [T21 

Direct and indirect dark matter detection rates are subject to uncertainties from both 
particle physics, through the microphysical properties of dark matter, and astrophysics. 
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FIG. 11: Flowchart illustrating the possible implications of comparing i^hep, the predicted dark 
matter thermal relic density determined from high energy physics, and r^cosmoj the actual dark 
matter relic density determined by cosmological observations. 




Annihilation 




Scattering 



FIG. 12: Efficient annihilation, corresponding to large indirect detection rates, is related to efficient 
scattering, corresponding to large direct detection rates. 



through the spatial and velocity distributions of dark matter. If completed, the research 
program described in Sec. IIII Al to pin down the properties of WIMPs will effectively remove 
particle physics uncertainties. Dark matter search experiments then become probes of dark 
matter distributions. 

As an example, consider direct detection. Theoretical predictions of direct detection rates 
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FIG. 13: Theoretical predictions for the direct detection neutrahno-proton scattering cross section 
a, as a function of neutrahno mass m^^ = m^, for various mSUGRA models (dots) [s^, and the 
prediction of LCC2 (*). ILC studies will constrain the values of a and to be smaller than the 
extent of the -k plotting symbol [s^ . 



are given in Fig. UHl As is typically done in particle physics studies, a simple dark matter 
halo profile is assumed throughout this figure. The enormous variation in rates results 
from particle physics uncertainties alone. LHC and ILC studies will reduce this uncertainty 
drastically. For example, for LCC2, the dark matter mass will be determined to a GeV 
at the ILC, and the cross section for neutralino-proton scattering will be determined to 
Act /a ^ 10% [s^. This constraint is shown in Fig. where the uncertainties are smaller 
than the extent of the * plotting symbol. 

Once collider constraints effectively remove microphysical uncertainties, the direct detec- 
tion rates give us information about the local dark matter density and velocity profile. In 
a similar way, indirect detection rates will provide additional complementary information. 
For example, experiments such as HESS and GLAST may detect photons from dark matter 
annihilation in the galactic center. Such rates are sensitive to the halo profile at the galactic 
center, a quantity of great interest at present. The synergy between collider experiments 
and these dark matter experiments will constrain the phase space distribution of WIMP 
dark matter in the Universe. Together with A^-body simulations, semi-analytical analyses of 
galaxy formation, and cosmology observations, these data will have important implications 
for the formation and evolution of structure. 
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FIG. 14: In superWIMP scenarios, a WIMP freezes out as usual, but then decays to a super WIMP, 
a sup erweakly-inter acting particle that forms dark matter. 



IV. SUPERWIMPS 



In superWIMP scenarios [3J, |33[, a WIMP freezes out as usual, but then decays to a 
stable dark matter particle that interacts superweakly, as shown in Fig. El The prototypical 
example of a superWIMP is a weak-scale gravitino produced non-thermally in the late 
decays of a weakly- interacting next-to- lightest sup ersy mmetric particle (NLSP), such as a 
neutralino, charged slepton, or sneutrino js^, 3^, 5y, 57, H, 5^, 60, 61 1. Additional examples 
include axinos and quintessinos |^ in supersymmetry, Kaluza-Klein graviton and 

axion states in models with universal extra dimensions j6^, and stable particles in models 
that simultaneously address the problem of baryon asymmetry 6^. SuperWIMPs have 
all of the virtues of WIMPs. They exist in the same well-motivated frameworks and are 
stable for the same reasons. In addition, in many cases the WIMP and superWIMP masses 
have the same origin. In these cases, the decaying WIMP and superWIMP naturally have 
comparable masses, and superWIMPs also are automatically produced with relic densities 
of the desired order of magnitude. 

As noted above, superWIMPs exist in many different contexts. We concentrate here on 
the case of gravitino superWIMPs. In the simplest supersymmetric models, supersymme- 
try is transmitted to standard model superpartners through gravitational interactions, and 
supersymmetry is broken at a high scale. The mass of the gravitino G is 



(11) 
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and the masses of standard model superpartners are 



m ~ — , (12) 

where M, = (87rG^)-i/2 ^ 2.4 x 10^^ GeV is the reduced Planck scale and F ~ (10^^ GeV)^ 
is the supersymmetry breaking scale squared. The precise ordering of masses depends on un- 
known, presumably C(l), constants in Eq. (fT^ . As discussed in Sec. IIIH most supergravity 
studies assume that the LSP is a standard model superpartner, such as the neutralino. 

The gravitino may be the LSP, however. In supergravity with high-scale supersymmetry 
breaking, the gravitino has weak scale mass M„eak ~ 100 GeV and couplings suppressed by 
M*. The gravitino's extremely weak interactions imply that it is irrelevant during thermal 
freeze out. The NLSP therefore freezes out as usual, and if the NLSP is a slepton, sneutrino, 
or neutralino, its thermal relic density is again I^nlsp ~ 0.1. However, eventually the NLSP 
decays to its standard model partner and the gravitino. The resulting gravitino relic density 
is ^ _ 

~ ^^NLSP • (13) 

"^NLSP 

In supergravity, where ~ mNLSP^ the gravitino therefore inherits a relic density of the 
right order to be much or all of non-baryonic dark matter. The superWIMP gravitino 
scenario preserves the prime virtue of WIMPs, namely that they give the desired amount of 
dark matter without relying on the introduction of new, fine-tuned energy scales. 

The superWIMP scenario differs markedly from other gravitino dark matter scenarios 0, 
1^ . 15, 16, 1^, EM Ilk, 20., 66, 67, 68, 69]. In the earliest proposals, gravitinos were produced 
thermally at temperatures T ~ Mpi, with Qq ~ 0.1 obtained by requiring ~ keV. Such 
scenarios are disfavored now by the expectation of an intervening era of inflation, which 
would dilute such a primordial population. After inflation, however, gravitinos may be 
produced in an era of reheating. In this case, Qq ~ 0.1 is obtained for reheat temperatures 
^RH ~ 10^° GeV. In contrast to the superWIMP scenario, it is not clear that gravitino 
production during reheating has testable consequences, other than the existence of cold 
dark matter itself. In addition, the reheating scenario requires the introduction of a new 
scale, in contrast to the superWIMP production mechanism, where the relic density is a 
function of the Fermi and Planck scales only and is naturally in the desired range. It is 
important to note, however, that the reheating and superWIMP production mechanisms 
are not mutually exclusive. The current gravitino relic population may have components 
from both production mechanisms, resulting in a very simple scenario in which dark matter 
is composed of two populations of particles with different histories and effects on the early 
Universe |8l| . 

Because superWIMP gravitinos interact only gravitationally, with couplings suppressed 
by M*, they are impossible to detect in conventional direct and indirect dark matter search 
experiments. At the same time, the extraordinarily weak couplings of superWIMPs imply 
other testable signals. The NLSP is a weak-scale particle decaying gravitationally and so 
has a natural lifetime of 

^ ^ 10^ - 10« s . (14) 

weak 

This decay time, outlandishly long by particle physics standards, implies testable cosmolog- 
ical signals, as well as novel signatures at colliders. 
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A. Cosmology 



The most sensitive probes of late decays with hfetimes in the range given in Eq. ()14|) are 
from Big Bang nucleosynthesis (BBN) and the Planckian spectrum of the cosmic microwave 
background (CMB). The impact of late decays to gravitinos on BBN and the CMB are 
determined by only two parameters: the lifetime of NLSP decays and the energy released 
in these decays. The energy released is quickly thermalized, and so the cosmological signals 
are insensitive to the details of the energy spectrum and are determined essentially only by 
the total energy released. 

The width for the decay of a slepton to a gravitino is 



rrir 



487rM2 m% 

* G 



1 

2 



(15) 



assuming the lepton mass is negligible. (Similar expressions hold for the decays of neutralino 
NLSPs.) This decay width depends on only the slepton mass, the gravitino mass, and the 
Planck mass. In many supersymmetric decays, dynamics brings a dependence on many su- 
persymmetry parameters. In contrast, as decays to the gravitino are gravitational, dynamics 
is determined by masses, and so no additional parameters enter. In particular, there is no 
dependence on left-right mixing or flavor mixing in the slepton sector. For mQ/mj ^ 1, the 
slepton decay lifetime is 



t{I^IG) ~3.6x 10^ 



100 GeV' 



rrii — m. 



G 



4 , 



TeV 



(16) 



This expression is valid only when the gravitino and slepton are nearly degenerate, but it is 
a useful guide and verifies the rough estimate of Eq. (fT^ . 
The energy release is conveniently expressed in terms of 



^EM = CEM-BeM^NLSP 



(17) 



for electromagnetic energy, with a similar expression for hadronic energy. Here esM is the 
initial EM energy released in NLSP decay, and -Bem is the branching fraction of NLSP decay 
into EM components. Inlsp = 't-nlsp/'^7 NLSP number density just before NLSP 

decay, normalized to the background photon number density = 2({3)T^/7i'^. It can be 
expressed in terms of the superWIMP abundance: 



Ynlsp ^ 3.0 X 10 



-12 



TeV 



0.23 



(18) 



Once an NLSP candidate is specified, and assuming superWIMPs make up all of the dark 
matter, with = Q^m = 0.23, the early Universe signals are completely determined by 
only two parameters: rriQ and ^nlsp- 



1 . BBN Electromagnetic Constraints 



BBN predicts primordial light element abundances in terms of one free parameter, the 

7- 



baryon-to-photon ratio rj = riB/n^- In the past, the fact that the observed D, "^He, ^He, and 
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^Li abundances could be accommodated by a single choice of t] was a well-known triumph 
of standard Big Bang cosmology. 

More recently, BBN baryometry has been supplemented by CMB data, which alone yields 
^710 = Tj/lO"^^ = 6.1 ± 0.4 This value agrees precisely with the value of rj determined 
by D, considered by many to be the most reliable BBN baryometer. However, it highlights 
slight inconsistencies in the BBN data. Most striking is the case of ^Li. For rjiQ = 6.0 ± 0.5, 
the value favored by the combined D and CMB observations, the standard BBN prediction 
is FT 



^Li/H = 4.71°:^ X 10-^° (19) 
at 95% CL. This contrasts with observations. Three independent studies find 



^Li/H = 1.5^°:^ X 10-^° (95% CL) |n| (20) 
^Li/H = 1.72;°;| X 10'^° (la + sys) ^ (21) 
^Li/H = 1.23;°:^^ X 10-^° (stat + sys, 95% CL) ^ , (22) 

where depletion effects have been estimated and included in the last value. Within the 
published uncertainties, the observations are consistent with each other but inconsistent 
with the theoretical prediction of Eq. ()19p. with central values lower than predicted by a 
factor of 3 to 4. ^Li may be depleted from its primordial value by astrophysical effects, 
for example, by rotational mixing in stars that brings Lithium to the core where it may 
be burned |7^l75|, but it is controversial whether this effect is large enough and consistent 
with the relatively small scatter of observations to reconcile observations with the BBN 
prediction j73l |. 

We now consider the effects of NLSP decays to gravitinos. For weakly-interacting NLSPs, 
that is, sleptons, sneutrinos, and neutralinos, the energy released is dominantly deposited in 
electromagnetic cascades. For the decay times of Eq. p4|) . mesons decay before they interact 
hadronically. The impact of EM energy on the lig ht element abundances has been studied 
m Refs. FaliiHlTilicl. The results of Ref. l79| are given in Fig. ^1 The shaded regions 



are excluded because they distort the light element abundances too much. The predictions 
of the superWIMP scenario for a stau NLSP with rriQ and mNLSP varying over weak scale 
parameters are given in Fig. by the grid. 

We find that the BBN constraint excludes some weak scale parameters. However, much of 
the weak scale parameter space remains viable. Note also that, given the ^Li discrepancy, the 
best fit is not achieved at ^em = 0, but rather for r ~ 3 x 10^ s and ^em ~ 10~^ GeV, where 
''Li is destroyed by late decays without changing the other relic abundances. This point is 
marked by the circle in Fig. ^1 The energy release predicted in the superWIMP scenario 
naturally includes this region. The ''Li anomaly is naturally resolved in the superWIMP 
scenario by a stau NLSP with mNLSP ~ 700 GeV and ~ 500 GeV. 



2. BBN Hadronic Constraints 

Hadronic energy release is also constrained by BBN [s^, E^, 113, EH EH E3- In fact, 
constraints on hadronic energy release are so severe that even subdominant contributions to 
hadronic energy may provide stringent constraints. 
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T (sec) 



FIG. 15: Predicted and excluded regions of the (r, Cem) plane in the super WIMP dark matter 
scenario, where r is the hfetime for / — > IG, and Cem is the normalized electromagnetic energy 
release. The grid gives predicted values for rriQ = 100 GeV — 3 TeV (top to bottom) and Am = 
mj — = 600 GeV — 100 GeV (left to right), assuming = 0.23. BBN constraints exclude 
the shaded regions; the circle indicates the best fit region where '''Li is reduced to observed levels 
without upsetting other light element abundances. Contours of CMB ^ distortions indicate the 
current bound (/i < 0.9 x 10~^) and the expected future sensitivity of Diffuse Microwave Emission 
Survey (DIMES) (/x ~ IQ-^). From Ref. Q. 



Slepton and sneutrino decays contribute to hadronic energy through the higher order 

processes 

/ ^ IZG , uWG 

V vZG , IWG , (23) 

when the Z or W decays hadronically. These three-body decays may be kinematically 
suppressed when mi - — uIq < m]y,mz, but even in this case, four-body decays, such as 

I — * l'~f*G —>■ IqqG, contribute to hadronic cascades and may be important. The branching 



fractions for these decays have been calculated in Refs. [57|, |58[ . The end result is that these 
constraints are stringent and important, as they exclude regions of parameter space that 
would otherwise be allowed. At the same time, much of the parameter space in the case of 
slepton and sneutrino NLSPs remains viable. For details, see Refs. js^, 13 • 

In contrast to the case of slepton and sneutrino NLSPs, the neutralino NLSP possibility is 
very severely constrained by bounds on hadronic energy release. This is because neutralinos 
contribute to hadronic energy even through two-body decays 

X ^ ZG, hG , (24) 

followed hy Z,h ^ qq. The resulting hadronic cascades destroy BBN successes, and exclude 
this scenario unless such decays are highly suppressed. Kinematic suppression is not viable, 
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however — if — < mz, the decay x iG takes place so late that it violates bounds 
on EM cascades. Neutralino NLSPs are therefore highly disfavored ; they 

are allowed only when the two-body decays to Z and h bosons are suppressed dynamically, 
as when the neutralino is photino-like, a possibility that is not well-motivated by high energy 
frameworks. 

3. CMB Constraints 

The injection of electromagnetic energy may also distort the frequency dependence of the 
CMB black body radiation (8^^J. For the decay times of interest, with redshifts z ~ 10^ to 
10^, the resulting photons interact efficiently through 7e~ 7e~ and eX — >■ 6X7, where X 
is an ion, but photon number is conserved, since double Compton scattering 76" — > 776" is 
inefficient. The spectrum therefore relaxes to statistical but not thermodynamic equilibrium, 
resulting in a Bose-Einstein distribution function 

A(^) = _ 1 ' (25) 

with chemical potential 7^ 0. 

In Fig. El we show contours of chemical potential /x, as determined by updating the 
analysis of Ref . 8S] . (For a more recent analysis and its implications for superWIMPs, see 
Ref. 89].) The current bound is < 9 x 10~^ |91i]. We see that, although there are 
at present no indications of deviations from black body, current limits are already sensitive 
to the superWIMP scenario, and are even beginning to probe regions favored by the BBN 
considerations described above. In the future, the Absolute Radiometer for Cosmology, 
Astrophysics, and Diffuse Emission (ARCADE) and Diffuse Microwave Emission Survey 
(DIMES) experiments may improve sensitivities to /i !^ 2 x 10~^ js^]. ARCADE and DIMES 
will therefore probe further into superWIMP parameter space, and will effectively probe all 
of the favored region where the ^Li underabundance is explained by decays to superWIMPs. 

B. Colliders 

The study of superWIMP dark matter at colliders has elements in common with the 
study of WIMPs, but with key differences. It may also be divided into three (overlapping) 
stages: 

1. SuperWIMP Candidate Identification. Is there evidence for late decays to super- 
WIMPs from collider studies? 

2. SuperWIMP Relic Density Determination. What are the superWIMP candidates' 
predicted relic densities? Can they be significant components or all of dark matter? 
What are their masses, spins, and other quantum numbers? 

3. Mapping the SuperWIMP Universe. Combined with other astrophysical and cosmo- 
logical results, what can collider studies tell us about astrophysical questions, such as 
the distribution of dark matter in the Universe? 
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For Stage 1, collider evidence for superWIMPs may come in one of two forms. Collider 
experiments may find evidence for charged, long-lived particles. Given the stringent bounds 
on charged dark matter, such particles presumably decay, and their decay products may be 
superWIMPs. Alternatively, colliders may find seemingly stable WIMPs, but the WIMP 
relic density studies described in Sec. 1111 Al may favor a relic density that is too large, a 
conundrum that may be resolved by postulating that WIMPs decay. These two possibilities 
are not mutually exclusive. In fact, the discovery of charged long-lived particles with too- 
large predicted relic density is a distinct possibility and would provide strong motivation for 
superWIMP dark matter. 

In the following subsections, we will explore how well the LHC and ILC may advance 
Stages 2 and 3. 



C. Relic Density Determination 

SuperWIMPs are produced in the late decays of WIMPs. Their number density is there- 
fore identical to the WIMP number density at freeze out, and so, as noted in Eq. (fT^ . the 
superWIMP relic density is 

"sWIMP — "WIMP • {^0) 

'^WIMP 

To determine the superWIMP relic density, we must therefore determine the superWIMP's 
mass. This is not easy, since the WIMP lifetime may be very large, implying that super- 
WIMPs are typically produced long after the WIMPs have escaped collider detectors. 

For concreteness, consider the case of supersymmetry with a stau NLSP decaying to a 
gravitino superWIMP. (Recall that, if superWIMPs are produced in sufficient numbers to 
be much of the dark matter, neutralino NLSPs are heavily disfavored, as their late decays 



invariably violate constraints from BBN and the CMB |32l . l33l l57l l58l . l61|.) As discussed 
above, the stau's lifetime is outlandishly long by particle physics standards. This gravitino 
superWIMP scenario therefore implies that the signal of supersymmetry at colliders will be 
meta-stable sleptons with lifetimes of days to months. Given their large mass, some of these 
sleptons will be slow, and so will produce highly-ionizing tracks that should be spectacularly 
obvious at the LHC 0, H H |3 • 

At the same time, because some sleptons will be slowly moving and highly-ionizing, they 
may be trapped and studied H |99, lOOj. As an example, sleptons may be trapped 
in water tanks placed outside collider detectors. These water tanks may then be drained 
periodically to underground reservoirs where slepton decays may be observed in quiet en- 
vironments. This possibility has been studied in Ref. (98| and is illustrated in Fig. ^1 
Alternatively, sleptons may be trapp ed in the detectors themselves jo^, 9^, or may be 
stopped in the surrounding rock |lO0j. 

How many sleptons may be stopped in a reasonably sized trap? The answer is highly 
model-dependent. The results for one model with 219 GeV sleptons is shown in Fig. ^1 
At the LHC, hundreds of sleptons may be caught each year in a 10 kton trap, assuming a 
luminosity of 100 fb~^/yr. A 10 kton trap is not particularly bulky. The optimal shape is one 
that covers as much solid angle as possible and is only ~ 1 m thick j9^ . These LHC results 
may be improved significantly if long-lived NLSP sleptons are kinematically accessible at the 
ILC. For the identical case with 219 GeV sleptons, 0(1000) sleptons may be trapped each 
year in a 10 kton trap at the ILC, assuming 300 fb~^/yr. By considering the slightly more 
general possibility of placing lead or other dense material between the ILC detector and the 
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FIG. 16: Left: Configuration for slepton trapping in gravitino superWIMP scenarios. Right: The 
number of sleptons trapped per year at the ILC in 10 kton (sohd), 1 kton (dot-dashed), and 0.1 
kton (dashed) water traps. The total number of sleptons produced is also shown (upper dotted) 
along with the number of sleptons trapped in the ILC detector (lower dotted). The trap shape 
and placement have been optimized and a luminosity of 300 fb~^/yr is assumed. The underlying 
model is minimal supergravity with = 600 GeV, mo = 0, Aq = 0, tan/3 = 10, and ^ > 0. 

The LHC results for this model are as indicated I98II. 



slepton trap, a further enhancement of an order of magnitude may be possible, allowing up 
to (9(10"^) sleptons to be trapped per ILC year. These ILC results are made possible by the 
ability to tune the beam energy to produce slow NLSPs. The abihty to prepare initial states 
with well-known energies and the flexibility to tune this energy are well-known advantages 
of the ILC. Here, these features are exploited in a qualitatively new way to produce slow 
sleptons that are easily captured. 

If thousands of sleptons are trapped, the slepton lifetime may be determined to the few 
percent level simply by counting the number of slepton decays as a function of time. The 
slepton mass will be constrained by analysis of the collider event kinematics. A percent 
level measurement of the slepton lifetime given in Eq. (|15|) therefore implies a high precision 
measurement of the gravitino mass, and therefore a determination of the gravitino relic 
density through Eq. As with the case of WIMPs, consistency at the percent level 

with the observed dark matter relic density will provide strong evidence that dark matter is 
indeed composed of gravitino superWIMPs. 

SuperWIMP quantum numbers and couplings may also be determined through collider 
studies [101, 102], although, as indicated above, these will typically be determined after or 
at the same time as the relic density determination, in contrast to the case of WIMPs. For 
example, an alternative method to determine the gravitino mass is to measure the energy 
of slepton decay products. This provides a consistency check of the mass determination 
described above. Alternatively, these two methods, when combined, determine not only 
tuq, but also the Planck mass M^,. Given enough events, the gravitino spi n m ay also be 
constrained to be 3/2 through detailed analyses of angular distributions 
gravitino differentiated from other superWIMP candidates, such as the axino 



Olll . and the 
9^ . The spin 
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and couplings of the gravitino may therefore be determined, showing that the superWIMP 
is in fact the superpartner of the graviton and that nature is locally supersymmetric. 



D. Mapping the SuperWIMP Universe 

Collider studies of superWIMPs will have significant implications for the phase space 
distribution of dark matter. In fact, the discovery of superWIMPs may resolve current 
discrepancies and shed light on important and controversial issues in structure formation. 

In the standard cosmology, dark matter is assumed to be cold, as is the case with WIMP 
dark matter. Cold dark matter is remarkably successful in explaining the observed large scale 
structure down to length scales of ~ 1 Mpc. Despite its considerable virtues, however, cold 
dark matter appears to face difficulty in explaining the observed structure on length scales 
^ 1 Mpc. Numerical simulations assuming cold dark matter predict, for exam ple, o ver dense 



cores in galactic halos 103] and too many dwarf galaxies in the Local Group [104 [. 

T hese prob lems may be alleviated or resolved by superWIMP dark matter lOSl IiOGL 



1071 Il08l Il09| . SuperWIMPs are produced with relativistic velocities at late times, as we 
have seen. They therefore exhibit properties typically associated with warm dark matter, 
suppressing power on small scales and potentially resolving the problems of cold dark matter 
mentioned above. ^ The discovery of superWIMP dark matter and the determination of 
NLSP and superWIMP masses and other relevant parameters at colliders would therefore 
change fundamentally our understanding of how galaxies were formed and provide a new 
framework for understanding halo profiles and the distribution of dark matter. 

As discussed above, decays that produce superWIMPs also typically release electromag- 
netic and hadronic energy. This energy may modify the light element abundances pre- 
dicted by standard BBN [3^ |3^ [t^ [sO, or distort the black body spectrum of the 
CMB B H H H 113. Collider studies will be able to determine how much energy is 
released and at what time, providing still more information with important consequences 
for astrophysics and cosmology. 



V. CONCLUSIONS 

This is an exciting time at the boundary of particle physics and cosmology. While impor- 
tant microphysical questions related to electroweak symmetry breaking and fiavor remain, 
breakthroughs in cosmology have added a whole new layer of fundamental problems requir- 
ing particle physics answers. 

Many of the key problems revolve around the mysteries of dark matter. Although there 
are many viable proposals, we have considered particularly well-motivated candidates whose 
relic densities are tied to two known energy scales, the Fermi and Planck scales, and fall 
"coincidentally" in the desired range. There are two classes of dark matter candidates with 
this property: WIMPs and superWIMPs. In both of these scenarios, dark matter particles, 
and typically many others, are expected with masses at the Fermi scale Mp ~ 100 GeV. In 
the next few years, particle colliders will at last probe this scale. If there is new physics. 



^ Note that gravitino's produced thermally during reheating are cold, and do not differ from standard cold 
dark matter in their impact in structure formation. 
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experiments at the Tevatron, LHC, and ILC will likely discover it and study it in great 
detail. 

If dark matter is composed of WIMPs, the LHC, and particularly the proposed ILC 
will be able to determine the WIMP's properties, and may also pin down its thermal relic 
density. If these determinations match cosmological observations to high precision, then 
(and only then) we will be able to claim to have determined what dark matter is. Such an 
achievement will also lead, through synergy with direct and indirect dark matter searches 
and cosmological observations, to improved knowledge of dark matter distributions and the 
formation of large scale structure. 

If dark matter is composed of superWIMPs, the LHC and ILC will again play a crucial 
role. The signal of new physics will likely be long-lived charged particles, a spectacular sig- 
nature that may be evident even in the first year or two of LHC running. By trapping these 
metastable particles and watching them decay, the properties and relic density of super- 
WIMPs may also be determined, providing another opportunity to identify dark matter. In 
the superWIMP scenario, the identification of "warm" dark matter may also resolve current 
issues in structure formation, and will also have implications for Big Bang nucleosynthesis 
and the cosmic microwave background. 

If any of the ideas discussed here is realized in nature, the coming years of exploration 
will not only provide our first incisive look at the Fermi scale, but will also yield profound 
insights about the Universe, its contents, and its evolution. 
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